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Abstract 

The model-independent constraints on the Abehan Z' couphngs from 
the LEP data are applied to estimate the Z' production in experiments 
at the Tevatron and LHC. The Z' total and partial decay widths are 
analyzed. The results are compared with model-dependent predictions 
and present experimental data from the Tevatron. If we assume the l-2a 
hints from the LEP data to be a signal of the Abelian Z' boson, then the 
Tevatron data constrain the Z' mass between 400 GeV and 1.2 TeV. 



1 Introduction 

Searching for signals of new physics beyond the standard model (SM) is an 
essential part of experiments at modern colliders. New phenomena could be 
discovered through deviations of observed quantities from the predicted SM 
background. However, observables in experiments at hadron colliders can be 
calculated with significant theoretical uncertainties coming mainly from the par- 
ton distribution functions of initial states and complicated structure of hadronic 
final states. In this situation one can only hope to discover the most prominent 
signals in the most clear processes. This is the reason to pay attention to search- 
ing for resonances of new heavy particles decaying into lepton pairs. 

A neutral vector boson (Z' boson) is probably the most perspective inter- 
mediate state in scattering processes of quarks and leptons which could be dis- 
covered in the Tevatron and LHC experiments. At the parton level it appears 
in the annihilation channel, its mass is allowed to be of order 1 TeV by current 
experimental constraints, and it is a necessary component of popular grand uni- 
fication theories and other models with extended gauge sector (see [U [3J [3] for 
review). 

In general, the accurate description of Z' resonance requires to consider 
scattering amplitudes with intermediate virtual states. But if the resonance is 
a narrow one, then it can be described in a more simple way by a small number 
of convenient characteristics of the production and the decay of the particle. In 
this approach it is enough to set the Z' mass and width, the production cross- 
section, and the branching ratio into the final state. Supposing some numbers 
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for the Z' parameters in various estimates one could and, in principle, would 
take into account all the available experimental constraints on Z' from previous 
experiments. 

Of course, effects of Z' boson can be calculated in details for each specific 
model beyond the SM. Such estimates are widely presented in the literature 
[H [SJ . Some set of popular based models and left-right models is usually 
considered in this approach. However, probing the set we can still miss the 
actual Z' model. In this regard, it is useful to complement model-dependent Z' 
searching by some kind of model-independent analysis, i.e. the analysis covering 
a lot of models. Almost all of the usually considered models belong to the models 
with so-called Abelian Z' boson. In Ref. [71 18] we found the relations which 
hold in any model containing the Abelian Z' boson and satisfying the following 
conditions: 

• only one neutral vector boson exists at energy scale about 1-10 TeV, 

• the Z' boson can be phcnomcnologically described by the effective La- 
grangian [TJ [31 [3] at low energies, 

• the Z' boson and other possible heavy particles are decoupled at consid- 
ered energies, and the theory beyond the Z' decoupling scale is either one- 
or two-Higgs-doublet standard model (THDM), 

• the SM gauge group is a subgroup of possible extended gauge group of the 
underlying theory. So, the only origin of possible tree-level Z' interactions 
to the SM vector bosons is the Z-Z' mixing. 

These relations cover almost all of the usually considered set of models (see 
|9l [To] for details). They require the same Z' couplings to the left-handed 
fermion currents within any SM doublet and the universal absolute value of 
the Z' couplings to the axial-vector currents for all the massive SM fermions. 
The relations reduce significantly the number of unknown Z' parameters. This 
allows to constrain the parameters by existing experiments as well as to predict 
the quantities used in the analysis of the Tevatron and LHC experiments. 

Recently we summarized the information about Z' couplings to leptons and 
quarks which can be extracted from the LEP experiments [HI [IHl- The Z' 
coupling to axial-vector currents was constrained by both LEP I and LEP II 
fi^ ,T^T^ data. In different processes it shows hints at about la confidence 
level (CL) with the approximately same maximum-likelihood (ML) value. This 
value can be used in estimates of observables in the Tevatron and LHC ex- 
periments. As for the couplings to vector currents, the Z' coupling constant 
to electron can be constrained by the LEP II e+e~ data only. Although the 
backward scattering shows a signal at the 2a CL, the ML value is outside of 
the 95% CL interval calculated by the complete set of bins. In this situation 
we refrain from using that ML value in our estimates. Nevertheless, the vector 
coupling is constrained at 95% CL. The upper bound on the electron vector 
coupling agrees closely with the corresponding upper bound on the axial- vector 
coupling. This fact allows us to suppose the rest of vector couplings to be con- 
strained by the same value, since no evident signals were discovered in other 
scattering processes measured by the LEP collaborations. It is worth to note 
that all the conclusions derived from the LEP data are also valid if one considers 
the THDM as the low-energy theory instead of the usual minimal SM. 
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The main goal of the present paper is to obtain estimates for the Z' parame- 
ters used in searching for the narrow resonance by applying the LEP constraints 
on the Z' couplings. Both the minimal SM and the THDM will be considered 
as the low-energy theory. 

The paper is organized as follows. Sec. [2] contains a necessary information 
about Z' interactions at low energies, the relations between the Z' couplings 
and the limits on these couplings obtained from the LEP data. In Sec. |3]the Z' 
production cross-section at hadron colliders is estimated. The bounds on the 
total and partial decay widths arc presented in Sec. El In Sec. [S] we discuss 
the application of our results comparing them, in particular, with the Tevatron 
experimental data and model-dependent predictions for the Tevatron and LHC. 
The explicit Lagrangian used for the calculations is given in Appendix [Xj 

2 Theoretical and experimental constraints on 
the Z' couplings 

In this paper we discuss mainly the Z' couplings to the vector and axial- vector 
fermion currents described by the Lagrangian 

l^zji = Izjr Mz + l^afi) cos 00 + {vf + 7'«/) sin ^o] /, 
Cz,ff = iz;/V[(^'/+7'a/)cos0o-(^^^^+7'4^)sin0o]/, (1) 

where / is an arbitrary SM fermion state; a/ and Vf are the Z' couplings to the 
axial- vector and vector fermion currents; is the Z-Z' mixing angle; v^^, 
are the SM couplings of the Z-boson. Such a parametrization is suggested by a 
number of natural conditions. First of all, the Z' interactions of renormalizable 
types are to be dominant at low energies ~ mw ■ The non-renormalizable inter- 
actions generated at high energies due to radiation corrections are suppressed 
by the inverse heavy mass 1/mz' (or by other heavier scales l/A; <C \/mz') 
and therefore at low energies can be neglected. It is also assumed that the Z' 
is the only neutral vector boson with the mass ~ mz' ■ 

It is obvious that the Lagrangian ([1]) requires the Z' boson to enter the the- 
ory as a gauge field through covariant derivatives with a corresponding charge. 
This idea allows also to introduce Z' couplings to SM scalar and vector fields. 
Although the latter couplings are inessential in the analysis of the Z' production 
cross-section in fermion collisions, they contribute to the Z' width. We assume 
that the SU(2)i x U(1)f gauge group of the SM is a subgroup of the GUT 
group. In this case, a product of generators associated with the SM subgroup 
is a linear combination of these generators. As a consequence, all the structure 
constants connecting two SM gauge bosons with Z' have to be zero. Hence, the 
Z' interactions to the SM gauge fields at the tree level are possible due to a 
Z-Z' mixing only. 

We will consider both the SM and the THDM as the low-energy theory. The 
explicit Lagrangian describing Z' couplings to the SM fields can be found in 
Appendix [X\ 

The parameters a/, vj, and 9o must be fitted in experiments. In a par- 
ticular model, one has some specific values for them. In case when the model 
is unknown, these parameters remain potentially arbitrary numbers. In most 
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investigations they are usually considered as independent ones. However, this 
is not the case if one assumes that the underlying extended model is a renor- 
malizable one. In Refs. [3 [5] it was shown that these parameters are correlated 
as 

Vf -af = vf, - af, , af = T^jgY^, (2) 

where / and /* are the partners of the SU{2)l fermion doublet {I* = i^ijiy* = 
hlu — Id ^'^d = ^3/ is the third component of weak isospin, and gY^ 
determines the Z' interactions to the SM scalar fields (see Appendix for details) . 
The parameter gY^p defines also the Z-Z' mixing angle in ([T]). 

As it was discussed in [9l [10] , the relations ([2]) cover a popular class of models 
based on the Eg group (the so called LR, x-V" models) and other models, such as 
the Sequential SM [TT]. Thus, they describe correlations between Z' couplings 
for a wide set of models beyond the SM. That is the reason to call the relations 
model-independent ones. 

The couplings of the Abclian Z' to the axial-vector fermion current have 
a universal absolute value. The value is proportional to the Z' coupling to 
scalar fields. Then, the Z-Z' mixing angle can be also determined by the 
axial- vector coupling. 

At low energies the Z' couplings enter the cross-section together with the 
inverse Z' mass, so it is convenient to introduce the dimensionless couplings 

■mz _ mz , . 

\/ATrmz' VAttttiz' 

which are constrained by experiments. Since the axial-vector coupling is uni- 
versal, we will use the notation 

a = = fle- = -a« = -ajy. (4) 

Then the Z-Z' mixing is 

sin9w COS0W mz 
t/Q ~ —2a -= . (5) 

V"em mz' 

It also follows from ([2]) that for each fermion doublet only one vector coupling 
is independent: 

^'/d = ""U + 2a- (6) 
As a result, Z' couplings can be parameterized by seven independent constants 

a, Vu, Vc, Vt, Ve, Wp, Vr- 

Recently we obtained limits on Z' couplings from the LEP I and LEP II data 
[9j[T0]. We found some hints of Z' boson at l-2a CL. Namely, the constants a 
and Ve show non-zero ML values. The axial- vector coupling a can be constrained 
by the LEP I data (through the mixing angle) and by the LEP II e~^e~ — >■ 
yU+^^,T+r^ data. The corresponding ML values are very close to each other. 
This value 

= 1.3 X 10"^ (7) 

will be used in our estimates. The 95% CL interval was also obtained by the 
experimental data: 

< < 3.61 X 10"^. (8) 
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The electron vector coupling Vg can be constrained by the LEP II e+e~ — > 
e^e~ data. An evident non-zero ML value occurred in fits taking into account 
the backward scattering bins only. Those fits showed 2a signal of the Z' boson. 
On the other hand, that ML value was excluded at 95% CL by fits including all 
the bins. This instability is the reason to refrain from using the ML value of Ve 
in our estimates. The 95% CL interval on Ve will be taken into account only: 

4 X 10"^ <vl< 1.69 X 10"^ (9) 

Other Z' coupling constants cannot be severely constrained by existing 
data. Among them Wc, and play an important role in the process 
qq ^ Z' ^ which is most perspective to discover the Z' resonance. 

Taking into account that no evident signals of new physics were found by the 
LEP collaborations in the processes involving quarks, muons and tau-leptons, 
we constrain the values of Vu, Vc, Vt, Vfj_, and Vt by the widest interval from the 
95% CL intervals for a and Ve'. 

< ^^other / < 4 X lO-\ (10) 

The knowledge of possible values of the Z' couplings allows to estimate the 
Z' production cross-section at the LHC and Tevatron and the Z' decay width 
without specifying the model beyond the SM. 



3 Z' production cross-section 

In modern experiments Z' bosons are expected to be produced in proton- 
antiproton collisions pp — ^ Z' (Tevatron) or proton-proton collisions pp — > Z' 
(LHC). At the parton level both the processes arc described by the annihilation 
of a quark- antiquark pair, qq — > Z' (Fig. [ij. The Z' production cross-section 
is the result of integration of the partonic cross-section aqq^z' with the parton 
distribution functions: 

<^AB = y2 '^^1 / fq^A{Xq,Q'^)fq.B{Xq,Q'^) 

Xaqq^Z'{mz',XqkA,XqkB), (11) 

where A, B mark the interacting hadrons (p or p) with the four- momenta fc^, 
ks] fq.A is the parton distribution function for the parton q in the hadron A 
with the momentum fraction Xq (0 < a;^ < 1) at the energy scale Q^. In our case 
= m^'. We use the parton distribution functions provided by the MSTW 
PDF package [12]. 

The production cross-section is determined by quadratic combinations of the 
Z' couplings to quarks, 

(TAB = 5^0-52 + aw„crac„ + w^ctbJ + a.Vc<Tav, 

+v^(Jfj2 + avtCFavt + v1a^2 . (12) 

where relations (jH)-® are taken into account. The factors a depend on mz' , 
the process type (proton- proton or proton-antiproton collision), and the beam 
energy. The factors Oav^ , (Jy2 , aavt and (Tg2 are small compared to 0-52 , (Tav^ and 
(Tg2 and their contributions to the cross-section can be neglected. 
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1.96 TeV 




Figure 1: Z' production at 
the parton level. 




Figure 2: Z' production cross-section 
vs. raz' in collisions at \fS ~ 1.96 
TeV. The filled area corresponds to the 
95% CL estimate, and the hatched area 
is for the ML estimate. 



We take into account the 90% CL uncertainties of the parton distribution 
functions provided by the MSTW PDF package. Finally, the production cross- 
section reads: 

cr = a^(Ta2 -I- au„(Tac„ + ^u^cj ± AcrP'^^, 
AaPdf = a^^aPf + ai;„A^L'!+CA<f. (13) 

Due to the existence of the ML value for the axial- vector coupling we perform 
two different estimates for the production cross-section: 

• 95% CL estimate. In this scheme both the couplings a and Vu are varied 
in their 95% CL intervals (|5|), ([TU)) . Then the production cross-section 
lies inside of the interval between zero and some maximal value. The 
maximal value is reached when both the couplings a and w„ are of the 
same sign and take their maximal values: a = v3l31 x 10~^, = 0.02. 
The uncertainty from the parton distribution functions should be also 
added. This estimate leads to the widest interval of possible values of the 
production cross-section. 

• Maximum-likelihood estimate. In this approach the axial-vector coupling 
is substituted by its ML value a = Vl.3 x 10~^. The vector coupling w„ 
is varied in its 95% CL interval. If one chooses the positive value of the 
axial-vector coupling, then the minimal value of the cross-section corre- 
sponds to Vu — —0.02 whereas the maximal value is reached at u„ ~ 0.02. 
The obtained interval should be also enlarged by ActP'^^ This estimate 
gives a more narrow interval for the production cross-section which can 
be considered as an 'optimistic' scenario to discover the Z' boson. 

The estimates for the Z' production cross-section in proton-antiproton col- 
lisions at the Tevatron and in proton-proton collisions at the LHC are shown in 
Figs. [Hand 131 respectively. In the LHC case the value is taken to be 7 TeV 
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Figure 3: Z' production cross-section vs. mz' in p-p collisions at \fS = 7 TeV 
and a/^ = 14 TeV. The filled area corresponds to the 95% CL estimate, and 
the hatched area is for the ML estimate. 




Figure 4: One-particle-irreducible correction to Z' Z' . 



and 14 TeV, corresponding to the current and expected energies. The Z' mass 
is chosen to be from 600 to 980 GeV for the Tevatron process and from 800 
to 2000 GeV for the LHC processes. At these masses it is possible to perform 
direct searches, and the boson production rate is not suppressed by the parton 
density effects. 



4 Z' width 

The Z' decay width V z' can be calculated by using the optical theorem: 

Yz: = J-^Gjml,) ^ 
mz' 

where G{p^) is the two-point one-particle- irreducible Green's function corre- 
sponding to the diagram in Fig. 01 We compute T z' at the one-loop level with 
the help of the FeynArts, FormCalc and LoopTools software [13j[T4]. The Feyn- 
man diagrams with internal Z' lines as well as the Passarino-Veltman integrals 
of type A give no contribution to the result, since they are real. The rest of 
diagrams correspond to different channels of Z' decay. As a result, we obtain 
also all the partial widths corresponding to Z' decays into two SM particles. 

All the Z' couplings to the SM scalar and vector bosons can be determined 
by the universal axial- vector constant ay and can be constrained. Then the 
partial widths corresponding to Z' decays into scalar and vector bosons are 
proportional to a?. As for the fermionic decays, the width can be written in the 
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form 

^z'^ff = (ij^aj + afVfTafVf + vjr^2 . (15) 
The factors Tg^2, TafVf and r^2 are proportional to mz'. Expressing eq. 

through the constants ([3]) one can see that the width is proportional to m^, and 
quadratic combinations of couphngs a, Vf. Thus it is convenient to introduce 
quantity 



mz' 

which is independent of niz' in our estimates. 

To calculate F numerically one has to choose values of the unknown masses 
of the SM scalar particles. If the minimal SM is considered as the low-energy 
theory, the only unknown mass is the Higgs boson mass m/j. The modern 
constraints on its value indicate that it is quite heavy, > 114 GeV. The 
contribution to the decay width from the scalar sector appears to be two or 
three orders of magnitude lower than the leading contribution from the fcrmionic 
decay channel. So the decay widths calculated at different values of are 
practically indistinguishable. In this regard, we present the results obtained for 
rrih = 125 GeV. 

When the THDM is considered, the scalar sector has six free parameters that 
can be expressed in terms of the masses nih, mn, TnAo, "^ff± and the mixing 
angles tana, tan/3 (see Appendix |X] for details). Because of the large number 
of physical scalar fields the estimates for the Z' width within the THDM can 
deviate from the results obtained in the case of the minimal SM. In order to 
obtain the most significant difference, we choose and Aq to be as light as it 
is allowed by the LEP constraints [15], namely 

mH± = 81 GeV, m^^ = 92 GeV. (17) 

The h and H masses are set to 

rrih ^mH = 125 GeV (18) 

just like in the SM case. The dependence of T on the mixing angles is negligibly 
weak. We take tan/3 = 2, which respects the LEP constraints. The tana value 
is set to 0.75. 

The decay width is estimated in two schemes which are similar to the case 
of the production cross-section: 

• 95% CL estimate. In this scheme the coupling constants a and Vf are 
varied in their 95% CL intervals ([TU)) . The minimal value of the width 
is calculated dX a = Vu = Vn.r = 0, We = ±a/0.4 x 10~^. The maximal 



value is realized when all the couplings are at their maximal absolute 
values, a and Vu,c,t are of the same sign, while Ve.fi.r have the opposite 
sign with respect to a: a ~ ±\/3.61 x 10~^, Vu,c,t = ±0.02, v^^r = ^0.02, 



T%/r69 X 10-2. 



Maximum-likelihood estimate. We set a = yOAS x 10"^ and vary w/ in 
their 95% CL intervals. We choose the positive value of a, so the min- 
imum value of the width corresponds to Ve = VOA x 10"^ and Vf = 
—OfTafUf /'2Ty2 (/ = fi,T,u, c,t). The maximum value is reached at 

Vu,c,t = 0.02, V^^r = -0.02, Ve = ~VTM X lO^^. 
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The Z' width ([TS]) is plotted in Fig. [S]as the function of Ue- The minimal 
SM and the THDM lead to slightly different bounds depicted in Fig. H 



f, GeV SM f, GeV THDM 




-0.01 0.01 _ -0.01 0.01 . 



Figure 5: The Z' width (HH) versus for the SM and THDM cases. The filled 
areas represent the 95% CL estimate, whereas the hatched areas represent the 
ML estimate. The inner vertical dot-dashed lines stand for the minimum 95% 
CL value of from the special one-parameter fit of the LEP II data, the outer 
ones depict the maximum 95% CL value of Ve from the general two-parameter 
fit of the LEP II data. 



f, GeV 95% CL f, GeV ML 




Figure 6: The F estimates versus Ve for the SM and THDM cases. The filled 
areas represent the estimate for the THDM case, and the hatched areas represent 
the estimate for the SM case. The meaning of the vertical dot-dashed lines is 
the same as in Fig. [S] 

Since we chose the positive ML value of the axial- vector coupling a. we obtain 
asymmetric domain in the parameter space within the ML estimate as it is seen 
in Figs. [5] and ini This asymmetry arises from the term avj^a^^ in (|15p. Of 
course, the sign of a is not constrained by the experimental data, so the sign of 
the vector coupling should be considered as the relative sign with respect to the 
axial- vector coupling. For the electron vector coupling the 2ct hint was observed 
[TU] . This allows to exclude the area near = shown in the figures. 

Consider an example of usage of the obtained estimates. Let us assume that 
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the Z' mass is of order 1-2 TeV, for instance mz' =1.5 TeV, so Z' production 
rate in the LHC and Tevatron processes is non-negligible and the direct searches 
are possible. The ML value F sa 50 GeV leads to the total decay width V z' = 169 
GeV. Thus we can expect the Z' resonance compatible with the narrow width 
approximation (NWA), T^,/m^, = 0.013 ^ 1. However, one has to keep in 
mind that Tz' ~ mz' is not excluded at the 95% CL. The extremely narrow 
resonances with F^/ w 1 GeV are also not excluded. 

It is also useful to estimate the partial decay widths of the Z' boson. In this 
analysis we take the ML value of the axial- vector coupling a — vOJS x 10~^ 
and vary other couplings in their 95% CL intervals. The results are presented 
as the plots in which a partial width is depicted versus the total width. In this 
way the branching ratios can be easily obtained. 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 20 40 60 80 100 

r^V"' ^liadr. GeV 

(a) (b) 

Figure 7: (a) The ML estimates on the F versus Fe+e- and Fp+^- plane. The 
filled area is for dimuon channel, and the hatched area represents the diclcctron 
channel, (b) The ML domain on the F-Fhadr plane. 



The partial decay widths for the electron-positron, muon-antimuon. and 
quark-antiquark channels are shown in Fig. [T] On these plots, the difference 
between the SM and the THDM case is negligible. As it is seen, the branching 
ratio for the electron-positron decay channel can be expected in the wide interval 

0.004 < ^ e+e") < 0.21. (19) 

Here, the minimal value corresponds to = 0, whereas the maximal value 
is reached at Ve ~ — vL69 x 10~^. The significant difference between the 
estimates for Fe+e- and F^+p- is caused by the fact that the Z' vector coupling 
to electron is much better constrained by the LEP II data that the muon one. 
The decay into quark-antiquark pairs can be the dominant decay channel. The 
corresponding probability can amount to 98%. 

Considering the Z' partial widths, one can find a significant distinction be- 
tween the SM and THDM in the scalar sector. Since a is the only Z' coupling 
entering the scalar and vector contributions to Tz' , there is the ML value of the 
partial decay width into two SM bosons (vectors or scalars). In the SM case, 
fbosons = 0.27 GeV. In the THDM case, f bosons = 0.53 GeV. The corresponding 
branching ratios arc less than 2.5%. 
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5 Discussion 



The recent experiments at the LEP gave some hints of the Abehan Z' boson. 
Although these hints correspond to 68-95% CL, they can be used as a beacon 
showing the most optimistic scenario to find Z' boson with a mass near 1 TeV. 
It is interesting to speculate about the question how can those hints look like 
at Tevatron and LHC experiments. Taking the LEP ML value of the axial- 
vector coupling we can give predictions under the assumption that a signal 
of the Abelian Z' boson has been probably observed in the LEP data. This 
estimate, called ML scheme, represents the most bold expectations concerning 
the Abelian Z' boson. Of course, such predictions do not exclude Z' boson with 
weaker axial-vector couplings. 

On the other hand the 95% CL bounds on possible Z' couplings to the SM 
particles arc left behind the LEP experiments. Taking these bounds for all the 
Z' couplings we can exclude some values of the observables at hadron colliders. 
In this scheme the values outside of the predicted intervals are forbidden for 
the Abelian Z' boson. Being measured in experiments, such values have to be 
interpreted as a signal of new physics which is something else than the Z' boson. 
For example, considering the Z' width, we can expect T z' x (1 TeV/m^')"^ ~ 
10 — 150 GeV from the ML estimate, and we can think about the NWA for 
mz' < 2 TcV. On the other hand, only extremely narrow resonances, V z' x 
(1 TeV/m^/)'^ < 1 GeV, and extremely wide resonances, T z' x (1 TeV/m^')"^ > 
500 GeV, can be surely excluded at the 95% CL. Thus, waiting for a narrow 
Z' resonance at hadron colliders we have to keep in mind that a more rich Z' 
phenomenology is still allowed by existing data. 

Now let us present the ML estimate for the Drell-Yan cross-section for the 
Tevatron and LHC experiments. As it was mentioned, in this case the NWA can 
be applied and the Z' contribution to the cross-section of the (pp) II process 
reads a(j>'p {pp) — >■ Z') x BR[Z' — > iT) where the branching ratio can be extracted 
from the total and partial Z' decay widths. The experimental bounds on the Z' 
contribution to the Drell-Yan process at the Tevatron are available in P^[T71[T5] 
together with the predictions from popular Z' models. The comparison between 
those results and our ML estimate for cr{pp — > Z' — > e+e^, yu+/i^) is presented 
in Figs. [H We can conclude from both the DO and CDF limits that the Z' 
hints from the LEP data can be the Abelian Z' boson with the mass between 
400 GeV and 1.2 TeV. Our model-independent results cover all the popular Z' 
models. We can also conclude that the model- independent lower bound on the 
Z' mass is still about 400 GeV whereas the popular models give the lower bound 
of order 800 - 900 GeV. 

It is straightforward to carry out similar calculations for pp ^ Z' ^ 
processes at the LHC. The ML domains are presented in Figs. IHl The cross- 
section values are plotted for the Z' mass up to 2 TeV. For higher masses the 
validity of the NWA is not guaranteed even for the ML estimate. Let us compare 
the results to the ones presented in [T9|. In Fig. 3 of Ref. [T9| the number of 
pp ^ Z' ^ I'^l' events for 100 fb~^ of integrated luminosity at ~ 14 TeV 
versus mz' is plotted. The ML number of pp ^ Z' ^ e~^e~ or n~ events 
for this luminosity can be obtained by multiplying the cross-section values in 
the left plot in Fig. H] by 10^. It can be seen that all the model-dependent 
predictions from Ref. |19j are covered by the e"*"e~ ML domain. 

In Table 2 of Ref. [19] the mo del- dependent estimates for a{pp — >■ Z' — > 
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Figure 8: The comparison between the Tevatron results and the ML estimates 
of the Z' production in the Drell-Yan process at y/S = 1.96 TeV. In all plots 
the filled areas represent the ML estimates. The experimentally obtained upper 
limits on the Z' contribution are taken from [TBI [H HH): the expected and 
observed 95% CL upper limits are depicted by the dashed lines and line charts, 
respectively, and the hatched areas are the 2a standard deviation bands for 
the expected values. The predictions from the popular models [ini HZl US] are 
plotted as sohd red hues, the corresponding models are Z[, Z!,^^, Z^, Z'^, Z!^, 
Z'^ and SSM Z' from the left to the right. 

Z+Z^) X Tz' are presented, mz' is set to 1.5 TeV. The ML estimate for this 
observable is easy to calculate using Figs. [3] and [7] (a) as a {pp ^ Z') X f ' X 
{mz'/l TeV)3. We obtain 94±92 pb-GeV and 210.7±210.1 pb-GeV for e+e" 
and ^J.^ fJ-^ decay channels, respectively. One can see that the predictions for the 
Z'^ and models (487 ± 5 fb • GeV and 630 ± 20 fb • GcV) lie outside the ML 
interval for the dielectron channel case, and the Z'^ prediction is not covered 
by the dimuon channel estimate. This is because mz' — 1.5 TeV appears to be 
quite heavy to provide exact value of the axial-vector coupling from the LEP 
data as it is assumed in the ML scheme. Of course, the model-dependent results 
are covered by the 95% CL intervals and cannot be excluded by the LEP data. 

The model-independent relations for the Z' couplings give a good possibility 
to reduce the number of unknown Z' parameters. As a consequence, the Z' 
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a{pp^Z)xBR(Z^tl'), fb 7TeV a{pp^Z)xBR[Z^fr), pb 14 TeV 




Figure 9: The ML domain for a{pp — > Z') x BR{Z' e+e ) (hatched area) 
and a{pp Z') x BR{Z' yr) (fihed area) oX-JS^l TeV and 14 TeV. 

width and the production cross-sections of the processes at modern hadron 
coUiders can be estimated using the constraints on the Z' couphngs obtained 
from previous experiments at LEP. A combined analysis of the LEP, Tevatron 
and LHC data seems to be possible. 

Our new model-independent results are complementary to the usual model- 
dependent schemes. The predictions of all the popular Z' models agree with 
our model-independent bounds. 

Finally the Z' hints observed in the LEP data can be still hidden as the 
resonance in the Tevatron experiments. We can expect this Z' boson with the 
mass between 400 GeV and 1.2 TeV. 



A Lagrangian 

In this section we adduce the scalar, fermion, Yukawa and gauge sectors of the 
Lagrangian that is used for the calculations. 

Let 0i (i = 1, 2) be two complex scalar doublets: 

+ + h + ic, \ 

^/2 I' ('"^ 

where Vi marks corresponding vacuum expectation values, are complex fields, 
and 6i, Ci are real fields. By diagonalizing the quadratic terms of the scalar po- 
tential F(0i,02) one obtains the mass eigenstates: two neutral CP-even scalar 
particles, R and /i, the neutral CP-odd scalar particle, Ao, the Goldstone boson 
partner of the Z boson, X3, the charged Higgs field, if^, and the Goldstone field 
associated with the boson, yj^: 

a\ = cos /3 — iJ+ sin /3, a\ = _ff cos /3 + sin ^, 
ci = X3 cos/3 - Ao sin/3, C2 = Ao cos/3 + xa sin/3, 
6i = i7 cos a — ft. sin a, 62 = /icosa + iJ sina, (21) 



where 



tan/3=— , (22) 
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and the angle a is determined by the exphcit form of the potential V{(i>\ ,4'2)- For 
instance, the CP-conserving potential, which has only CP- invariant minima, 
can be used [Ml HI]: 

2 

i=l 

+X4{lm[4<f>2]f + X5i4<f>i)i4<f>2)- (23) 



It is consistent with the absence of the tree-level flavor-changing neutral currents 
(FCNC's) in the fermion sector. The corresponding value of a is [55] 

tan2a = - 7%(^^ + y . (24) 
A2W2 - •^I'^i 

The Z' coupling to the scalar doublets can be parametrized in a model 
independent way as follows [52]: 



i=l 



(25) 



where g, g', g are the charges associated with the SU(2)i, U(l)y, and the Z' 
gauge groups, respectively, ctq are the Pauli matrices, 

= ( '0 ' tl ) <^'=) 

is the generator corresponding to the gauge group of the Z' boson, and Y^. is 
the U(1)y hypercharge. The condition Y^^ = 1 guarantees that the vacuum is 
invariant with respect to the gauge group of photon. 

The vector bosons, A, Z , and Z' , are related with the symmetry eigenstates 
as follows: 

B ~^ A cos 6w — [Z cos 9^ — Z' sin 6'o) sin 0vKj 
W-i Asin^vv + (^cos^o — ^'sin0o)cos0VKj 

B Z sin 6*0 + cos 6*0, (27) 

where XwlQw = <?'/<? is the adopted in the SM value of the Weinberg angle, and 

.gm^ (^01, 2 cos^ /3 -t- Y^^^2 sin^ ;3) 

tan 00 = \ , — -— ^. (28) 

g cos By/ (m^/ — / cos-^ 9^ ) 

As is seen, the mixing angle is of order ~ TO^/m|,. That results in the 
corrections of order ~ rn^/rn^, to the interactions between the SM particles. 
To avoid the tree-level mixing of the Z boson and the physical scalar field ^0 
one has to impose the condition Y^i,2 = Y^2,2 = ^0,2- 

The effective low-energy Lagrangian of the fermion- vector interactions reads 
[201 [211 [22]: 

fL 

E {^^ - - jB^Ynj) fn, (29) 
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where the renormahzable type interactions are admitted and the summation 
over the all SM left-handed fermion doublets, /l = {{fu)L, {fd)L}, and the 
right-handed singlets, fa = (/«)/£, {fd)R, is understood. Qf denotes the charge 
of / in the positron charge units, 

tf. - f¥ (») 



It,/, 

and Yf^ equals to —1 for leptons and 1/3 for quarks. 

In the present paper we use the Z' couplings to the vector and axial- vector 
fermion currents defined as 

^/ = -9^^^ «/ = ■ (31) 

The Lagrangian leads to the interactions between the fcrmions and the Z 
and Z' mass eigenstates described by (P). 

Renormahzable interactions of fermions and scalars are described by the 
Yukawa Lagrangian. To avoid the existence of the tree-level FCNC's one has to 
ensure that at the diagonalization of the fermion mass matrix the diagonalization 
of the scalar- fermion couplings is automatically fulfilled. In this case the Yukawa 
Lagrangian, which respects the SU(2)i x U(l)y gauge group, can be written in 
the form: 



z 

fL i=l 

+ GU,^ [M'lifu)R + {fu)Rcl^ffL\ } , (32) 



where 0^ = i<J24'i the charge conjugated scalar doublet, and the Cabibbo- 
Kobayashi-Maskawa mixing is neglected. Then, the fermion masses are 

TO/ = '^^^ {Gf^i cos/3 -1- Gf^2 sin/3) . (33) 

As was shown by Glashow and Weinberg [35], the tree-level FCNC's medi- 
ated by Higgs bosons arc absent in case when all fermions of a given electric 
charge couple to no more than one Higgs doublet. This restriction leads to four 
different models, as discussed in Ref. [2S]. In what follows, we will use the most 
general parametrization ((32)) including the models mentioned as well as other 
possible variations of the Yukawa sector without the tree-level FCNC's. 

The gauge sector is taken to be 



r — —-F'^'^F —-F^^^F"- —-Ft^'^F f34~l 



In the present paper all calculations are carried out in the Feynman-'t Hooft 
gauge, the gauge-fixing functions are 
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G 



( 



1 




(35) 



Then, the gauge-fixing part of the Lagrangian reads 



C, 



'gauge fixing 




(36) 



The kinetic part of Faddeev-Popov sector is 



'ghost 



uz'id^ + £,m%,)uz', 



+ 



(37) 



^ is gauge- fixing parameter. For arbitrary ^ the gauge-boson propagator is 



The MSM parametrization can be obtained by putting tan/3, tana, /i2, 
•^2.3,4.5, ^02,1,2: Gd.2 and Gu.2 to zero and dropping the summations over i in 
dsn) and dSni). 

In Refs. [71 [S] the relations between Z' parameters were found from the 
requirement that the underlying extended model is a renormalizable one. They 
read 



%,i^%.2 = %, Ylj^Ylj*, Yrj ^YLj + 2T3fY^ (39) 



in case of Abelian Z' boson. Here / and /* are the partners of the SU{2)l 
fermion doublet {I* = i>i,v* = l,qu = Qd and = (/„), T3/ is the third compo- 
nent of weak isospin. These relations are used all over the present paper. They 
can be also rewritten in terms of vector couplings, axial- vector couplings, and 
the Z-Z' mixing angle. 
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